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We present improved measurements of CP-violation parameters in B° — > 4>K°, f]'K°, KgKgKg 
decays based on a sample of 535 x 10 6 BB pairs collected at the T(4S) resonance with the Belle 
detector at the KEKB energy- asymmetric e + e~ collider. We obtain sm2<j>1 a = +0.64±0.10(stat) ± 
0.04(syst) for B° -> n'K°, +0.50 ± 0.21(stat) ± 0.06(syst) for B° -» cj>K° , and +0.30 ± 0.32(stat) ± 
0.08(syst) for B° -> K° S K° S K% decays. We have observed CP violation in the B° -> n'K° decay 
with a significance of 5.6 standard deviations. We also perform an improved measurement of CP 
asymmetries in B° -> J/i>K° decays and obtain sin20i = +0.642 ± 0.031(stat) ± 0.017(syst). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 



Particles from physics beyond the standard model 
(SM) may contribute to B° meson decays mediated by 
flavor-changing b — y s transitions via additional quan- 
tum loop diagrams, and potentially induce large devia- 
tions from the SM expectation for time-dependent CP 
asymmetries Q. In the decay chain T(45) -> B°B° -> 
/cp/tag, where one of the B mesons decays at time tcp 
to a CP eigenstate fc p and the other decays at time i tag 
to a final state /t ag that distinguishes between B° and 
£>° , the decay rate has a time dependence Q given by 



V(At) 



-|At|/r B o 
4r B o 



Sf sin(ATOdAi) 



A f cos(Am d At) \. (1) 



Here Sf and Af are CP- violation parameters, t b q is the 
B° lifetime, Am^ is the mass difference between the two 
B° mass eigenstates, At — tcp — itag, and the &- flavor 
charge q = +1 (—1) when the tagging B meson is a 
B° (TP). In the SM, CP violation arises only from the 
irreducible Kobayashi-Maskawa phase Q in the weak- 
interaction quark- mixing matrix jjjj. The SM approx- 



imately predicts Sf = — £/sm20i|5( and Af = for 
both b — y ccs and b —y sqq transitions, where £/ = +1 
(—1) corresponds to CP-even (-odd) final states, in the 
leading order. Recent SM calculations [6( for the effec- 
tive sin2</>! values, sin.2(j>\ s , obtained from B° — > (j>K°, 
rj'K and KgKgKg agree with sin 20! , as measured in 
B° —y J/yjK° decays, at the level of 0.01. Thus compar- 
ison of measurements of Sf and Af between modes is an 
important test of the SM. 

Previous measurements of CP asymmetries in b — y sqq 
transitions by Belle and BaBar Q differed from the 
SM expectation, although the deviation was not statis- 
tically significant. BaBar has since updated their re- 
sults [9]. In this Letter we describe improved mea- 
surements of S f and A f in B° -» <f>K%, (t)K Q L , rfK% 
and KgKgKg decays using a data sample of 492 fb _1 
(535 x 10 6 BB pairs), which is nearly twice that used for 
our previous measurements. The analysis has also been 
improved by the addition of the following decay chains: 



B° -> <t>K% 



S) 



K%Kl; B° 



V'K° L ; 



7r u 7r . We also describe improved measurements 
of sin 20i and Af in B° — * J/ipKg and J/ipK® decays 
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using the same data sample; our previous measurement 
used a 140 fb _1 data sample [10|. These modes have 
the largest statistics coupled with the smallest theoreti- 
cal uncertainties and thus provide a firm reference point 
for the SM. 

At the KEKB e nerg y-asymmetric e + e~ (3.5 GeV on 
8.0 GeV) collider [IJ, the T(4S) is produced with a 
Lorentz boost of /3j — 0.425 nearly along the elec- 
tron beamline (z). Since the B° and B° mesons are 
approximately at rest in the T(45) center-of-mass sys- 
tem (cms), At can be determined from the displace- 
ment in z between the fcp and /t ag decay vertices: 
At ~ (z CP - z tag )/(/3 7 c) = Az/(/3 7 c). 

The Belle detector [12j is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array 
of aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to 
detect K® mesons and to identify muons (KLM). 

Charged tracks reconstructed with the CDC, except for 
tracks from K s — > 7r + n~ decays, are required to originate 
from the interaction point (IP). We distinguish charged 
kaons from pions based on a kaon (pion) likelihood Ck(k) 
derived from the TOF, ACC, and dE/dx measurements 
in the CDC. Photons are identified as isolated ECL clus- 
ters that are not matched to any charged track. Can- 
didate K^ L mesons are selected from ECL and/or KLM 
hit patterns that are consistent with the presence of a 
shower induced by a K^ meson. 

The intermediate meson states are reconstructed from 
the following decays: it — > 77, K s — > ir + ir~ (de- 
noted by K%~) or 7r°7r° (denoted by K s °), 



ir^ir~, T]' — > p°7 or T]ir + ir~ , <f> — > K + K~ 
and J/ij) — > £ + £~ (£ = /i, e). We use all combina- 
tions of the intermediate states with the exception of 
{rj — » -K + ir~ir°, rf — » pj,K s °} candidates for B Q — > 
{r]'Kg ,r]'Kl,J/ipKg} decays, respectively. The recon- 
struction and selection criteria for B° — > fcp decay can- 
didates are almost the same as in our previous measu- 
ment 0, E3]. The K° L and K s ° candidates for rj'K 
are reconstructed with the same method as used for 

K s K l K s deca y 



V 



77 or 



4>K° [l]|. We reconstruct the B° 



in the Kg 



or K7 



-K^~K S ° final states. 



In 



Kg K L decays are used for the B c 



addition, 4 

4>Kg~ sample. We identify candidate B° — » fcp de- 
cays without a K*l meson using the energy difference 



AE = El 



i 'cam 



and the beam-energy constrained 

is the 



where ££ eam 



mass M hc = y/{EL a J 2 ~ WbJ 
beam energy, and E* B and p* B are the energy and momen- 
tum, respectively, of the reconstructed B candidate, all 
measured in the cms frame. The signal candidates are 
selected by requiring M^ c G (5.27,5.29) GeV/c 2 and a 



mode-dependent AE window. Only Mb c is used to iden- 
tify the decay B° -> <j)K s followed by <f> -> K S K° L . Other 
candidate B° — > fcp decays with a K^ are selected by 
requiring p* B G (0.2,0.45) GeV/c for B° -> J/i)Kl can- 
didates and p* B G (0.2, 0.5) GeV/c for the others. 

The dominant background for the b — > sqq signal 
comes from continuum events e + e _ — > qq where q = 
u, d, s, c. To distinguish these topological^/ jet-like events 
from the spherical B decay signal events, we combine a 
set of variables that characterize the event topology 
into a signal (background) likelihood variable C s - lg (£bkg), 
and impose loose mode-dependent requirements on the 
likelihood ratio TZ s /h = C sig /(C Big + £bk g )- 

The contributions from BB events to the background 
for B° — > fcp candidates with a K^ are estimated with 
Monte Carlo (MC) simulated events. The small (~ 3%) 
BB combinatorial background in B° —> rj'K s (r]' —> p°j) 
is estimated using MC events. We reject K s K s Kg candi- 
dates if one of the K s pairs has an invariant mass within 
±2cr of the XcQ mass or ±lcr of the D° mass, where a 
is the K S K S mass resolution. The fraction of B° — » 
K+K~K S and B° 



fo(980)K s (/ (980) -» K+K~) 
events in the B° —> <f>K s sample is estimated to be 
2.75 ± 0.14% and zero within error, respectively, from 
the Dalitz plot for B -> K+ K~ K candidates [lj. 

The 6-flavor of the accompanying B meson is identified 
from inclusive properties of particles that are not associ- 
ated with the reconstructed B° — > fcp decay. The tag- 
ging information is represented by two parameters, the 
6-flavor charge q and r 15j . The parameter r is an event- 
by-event, MC-determined flavor-tagging dilution factor 
that ranges from r = for no flavor discrimination to 
r = 1 for unambiguous flavor assignment. For events 
with r > 0.1, the wrong tag fractions for six r inter- 
vals, Wi (I = 1,6), and their differences between B° and 
B° decays, Awi, are determined using semileptonic and 
hadronic b — > c decays 0, Gil ■ If r < 0.1, we set the 
wrong tag fraction to 0.5, and therefore the tagging in- 
formation is not provided. The total effective tagging 
efficiency is determined to be 0.29 ± 0.01. 

The vertex position for the fcp decay is reconstructed 
using charged tracks that have enough SVD hits [lj|. 
The /tag vertex is obtained with well-reconstructed tracks 
that are not assigned to fcp- A constraint on the 
interaction-region profile in the plane perpendicular to 
the beam axis is also used with the selected tracks. 

Figures [lja)-(m) show the distributions of recon- 
structed variables 1Z s /b, M^ c and p* B after flavor tag- 
ging and vertex reconstruction for B° — > r)'K°, 4>K° 
and K S K S K S candidates. The M^ c distribution for the 
B° — > J/i/jK s candidates and p* B distribution for the 
B° — » J/ipK® candidates are shown in Fig. [2j The sig- 
nal yield for each mode is obtained from an unbinned 
maximum-likelihood fit to these distributions; the AE 
distribution is also included in the fit for the modes with- 
out a K L meson. The signal shape for each decay mode 
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Mbc and p* B distributions for reconstructed 
/b, Mbc with 72. s /b < 0.5 and Mb c with 



72. s /b > 0.5 distributions for (a, b and c) B° — > rj'Kg, (d, e and 
f) B° -> ^JC§, and (g, h and i) 5° -» 72 s/b andp^ 

for (j and k) B° -» 7j'Jf£ and (1 and m) B° -> The solid 

curves and histograms show the fits to signal plus background 
distributions, and hatched areas show the background contri- 
butions. Background contributions are subtracted in figures 
(k) and (m). 
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FIG. 2: (a) Mbc distribution in the AE signal region for se- 
lected B° — > J/ipKg candidates and (b) p B distribution for 
selected B° — > J/ipK ( [ candidates. 



is determined from MC events. The background has two 
components: continuum, which is modeled using events 
outside the signal region, and BB background, which 
is modeled with MC events. The signal yields are de- 
termined to be 307 ± 21 for B a -> <f>K% 114 ± 17 for 
B° -> <f>Kl, 1421 ± 46 for B° 
B" 



rj'K 



£, 185 ±17 for B° 



r]'K a s , 454 ± 39 for 



K s K s K S' 7484 ±87 for 



errors are statistical only. 

We determine Sf and Af for each mode by perform- 
ing an unbinned maximum-likelihood fit to the observed 
At distribution. The probability density function (PDF) 
for the signal distribution, V B i s (At; Sf,Af, q, wi, Awi), is 
given by EclJT]) fixing Tgo and Ama at their world aver- 
age values [Uj and incorporating the effect of incorrect 
flavor assignment. The distribution is convolved with the 
proper-time interval resolution function R s i g (At), which 
takes into account the finite vertex resolution [l6|]. We 
determine the following likelihood for each event: 

Pi = (1 - /oi) E / [Pk(At')R k (AU - At')] d(At') 



folPol(Ati), 



(2) 



where k denotes signal, continuum and BB background 
components. In the B° — » KgKgKg and J/ipKg samples 
the BB component is negligibly small and not included 
in the fit. The fraction of each component fk depends on 
the r region and is calculated on an event-by-event basis 
as a function of the following variables: AE and Mbc 
for B° -> J/i)K%- p* B for B° -> J/ipKl; p* B and K s/h 
for B° -> rj'Kl and ^K° L ; M h 



and TZ s / h for B° 



K%Kl)K s ] AE, M hc and K s/h for the other modes. The 
PDF Vk(At) for background events is convolved with the 
resolution function for the background 0,0. The 
term P ol (At) is a broad Gaussian function that represents 
a small outlier component with a fraction of f \ [l6| . The 
only free parameters in the fits are Sf and Af, which are 
determined by maximizing the likelihood function L = 
Yli Pi(Ati;Sf, Af) where the product is over all events. 

Table U summarizes the fit results for sin 2<j)f and Af. 
These results are consistent with and supersede our pre- 
vious measurements 0, E3] • Fits to each individual mode 



TABLE I: Results of the fits to the At distributions. The first 
errors are statistical and the second errors are systematic. 



Mode 


sin 2<f>f 


Af 


<f>K u 


+0.50 ±0.21 ±0.06 


±0.07 ±0.15 ±0.05 


T)'K° 


±0.64 ±0.10 ±0.04 


-0.01 ±0.07 ±0.05 


K%K%K% 


±0.30 ±0.32 ±0.08 


±0.31 ±0.20 ±0.07 


J/ipK° 


±0.642 ±0.031 ±0.017 


±0.018 ±0.021 ±0.014 



B" -> J/i>K% and 6512 ± 123 for B° -> J/tpK° L , where 



with K° s and K° L yield (S v>K o, A ri , K o) = (±0.67 ±0.11, 
-0.03 ±0.07), (5,,^, A n , K o) = (-0.46 ±0.24, ±0.09 ± 
0.16), (5 0Kg) A^ K o L ) = (±0.50 ± 0.23, ±0.11 ± 0.16), 
(S*K°> A*Kl) = (-0.46 ±0.56, -0.15 ±0.38), {S J/ilK o, 
Aj/^ko) = (±0.643±0.038, -0.001±0.028) and (S J/ilK o , 
A J/4 , K l) = (-0.641±0.057, ±0.045±0.033), where errors 
are statistical only. We define the background-subtracted 
asymmetry in each At bin by (N+ - N-)/(N+ + N-), 
where N+ (N-) is the signal yield with q = ±1 (—1). Fig- 
ures EKa)- (d) show the At distributions and asymmetries 
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FIG. 3: Background subtracted At distributions and asym- 
metries for events with good tags (r > 0.5) for (a) B° — > rj'K , 
(b) B° -> cj,K°, (c) S° -» XgA-gTfg and (d) B° -> J^ff . In 
the asymmetry plots, solid curves show the fit results; dashed 
curves show the SM expectation from our B° — > J/tpK mea- 
surement. 



for good tag quality (r > 0.5) events. The sign of each 
At measurement for final states with a K\ is inverted in 
order to combine results with Kg and K® mesons. 

The dominant sources of systematic error for sin 2</>° ff 
in b — > sqq modes come from the uncertainties in the res- 
olution function for the signal (0.03 for the B° — > rj'K 
mode, 0.04 for the (f>K° mode, 0.05 for the B° -> 
if^iCjiCy mode) and in the background fraction (0.02, 
0.04, 0.06). The effect of / (980)iC background in the 
(j)K° mode (0.02) is estimated using the BES measure- 
ment of the /q(980) lineshape [l8| and is included in 
the background fraction systematic error. The domi- 
nant sources for At in b — > sqq modes are the effects 
of tag-side interference [lj| (0.02, 0.03, 0.04), the un- 
certainties in the background fraction (0.02, 0.03, 0.05), 
in the vertex reconstruction (0.02 for all modes) and in 
the resolution function (0.02, 0.01, 0.02). We study the 
possible correlations between TZ s /b, p* B and r PDFs used 
for c^K^ and r/'K^, which are neglected in the nominal 
result, and include their effect in the systematic uncer- 
tainties in the background fraction. Other contributions 
come from uncertainties in wrong tag fractions, the back- 
ground At distribution, t b o and Am ( i- A possible fit bias 
is examined by fitting a large number of MC events and 
is found to be small. 

The dominant sources of systematic errors for the 
B° — » J/ipK° mode are the uncertainties in the vertex 
reconstruction (0.012 for sin2</>i, 0.009 for Af), in the 
resolution function for the signal (0.006, 0.001), in the 
background fraction (0.006, 0.002), in the flavor tagging 



(0.004, 0.003), a possible fit bias (0.007, 0.004) and the 
effect of the tag-side interference (0.001, 0.009). Other 
contributions amount to less than 0.001. We add each 
contribution in quadrature to obtain the total system- 
atic uncertainty. 

For the B° — >• rj'K mode, we observe CP violation 
with a significance equivalent to 5.6 standard deviations 
for a Gaussian error, where the significance is calculated 
using the Feldman-Cousins frequentist approach [2(ij |. 
The results for B° -> rj'K , (f>K° and K^K^K° S decays 
are all consistent with the value of sin 2(f>i obtained from 
the decay B° — > J/ipK° within two standard deviations. 
No direct CP violation is observed in these decay modes. 
Further measurements with much larger data samples are 
required to search for new, beyond the SM, CP-violating 
phases in the b — > s transition. 
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